Heat conduction and diffusion of hard disks in a narrow channel 
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Using molecular dynamics we study heat conduction and diffusion of hard disks in one dimensional 
narrow channels. When collisions preserve momentum the heat conduction re diverges with the 
number of disks N as re ~ N a (a ~ 1/3). Such a behaviour is seen both when the ordering of disks 
is fixed ('pen-case' model), and when they can exchange their positions. Momentum conservation 
results also in sound-wave effects that enhance diffusive behaviour and on an intermediate time 
scale (that diverges in the thermodynamic limit) normal diffusion takes place even in the 'pen-case' 
model. When collisions do not preserve momentum, re remains finite and sound-wave effects are 
absent. 
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According to the Fourier's law of heat conductivity, 
when a small temperature difference is applied across a 
system, in the steady state the heat flux j satisfies the 
equation 
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where T is a local temperature and re is the heat conduc- 
tivity of a given material. Since the heat conductivity 
re is one of the most important transport coefficients a 
lot of efforts were devoted to its calculation. Of particu- 
lar interest is the case of low-dimensional systems, where 
re might diverge and thus Eq.([T]) would break down 
However, despite intensive efforts the nature of this diver- 
gence is not yet fully understood even in one-dimensional 
systems. Renormalization group calculations [2] show 
that for one-dimensional fluid-like systems re should di- 
verge with the number of particles N as re ~ N a with 
a = 1/3. Although earlier works suggested other values 
of a, recent simulations for hard-core particle systems 
agree with this prediction [3J]. In the other class of sys- 
tems, chains of nonlinear oscillators (Fermi-Pasta-Ulam 
systems), simulations suggest [J H[ a larger value of a 
(0.37-0.40) and that might be consistent with the pre- 
dictions of mode-coupling theory a = 2/5 [E @1- How- 
ever, recently it was suggested that for such chains of 
oscillators a also should be equal to 1/3 and numerically 
observed values of a were attributed to numerical diffi- 
culties [Tj. 

It would be desirable to relate the divergence of heat 
conductivity to other dynamic properties of a given sys- 
tem. For example, it has been shown that chaoticity 
plays an important role and ensures that re remains fi- 
nite Q . On the other hand, conservation of momentum 
is known to imply the divergence of re [1, 0] ■ Some at- 
tempts were also made to relate heat conductivity and 
diffusion. In particular, Li and Wang suggested Q that 
the exponent (3 describing the mean square displacement 
of diffusing particles 
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should be related with a through the equation 
a = 2 - 2/(3. 



(2) 
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Such an equation implies that the normal diffusion (j3 = 
1) leads to the normal (non-divergent) heat conductivity 
(a = 0). Moreover, superdiffusion (f3 > 1) and subd- 
iffusion (/3 < 1) correspond to divergent (a > 0) and 
vanishing (a < 0) heat conductivity,respectively [lfj. 
Although some numerical examples [ll| seem to confirm 
the relation ([3]), its derivation is based only on quali- 
tative arguments that neglect, for example, interactions 
between particles and so the suggestions that the rela- 
tion J3J) is of more general validity should be taken with 
care [12j. In another attempt, studying a class of non- 
interacting billiard heat channels, Denisov et al. [IH ob- 
tained a different relation between the exponents a and 
j3, namely 



a = (3-1. 



(4) 



The relation ((4]) was verified numerically for the energy 
diffusion in a one dimensional hard-core model [141 ] . How- 
ever, it was argued [l4[ that the Levy walk scenario, that 
the energy diffusion obeys in this model, might be due 
the absence of exponential instability and it is not clear 
whether this result can be extended to more realistic sys- 
tems. 

Establishing a firm relation between heat conductiv- 
ity and diffusion could be possibly very influential and 
shed some light also on other transport phenomena. In 
the present paper we examine a model of hard disks in 
a narrow channel. When a fraction of disks is immo- 
bile and thus collisions do not conserve momentum, the 
heat conductivity re remains finite and normal diffusion 
takes place. Simulations show that when there are no 
immobile disks and momentum is conserved, heat con- 
ductivity diverges with a close to 1/3. According to ^ 
or ^ it should imply a superdiffusion. Although in this 
case diffusion is enhanced by sound-wave effects, there 
are no indications of superdiffusivity. Our work suggests 
that in hard-disk systems in a narrow channel heat con- 
ductivity and diffusion might be related but in a more 
intricate way than relations ([3]) or Q would suggest. 

In our model N identical hard disks of radius r and 
unit mass are moving in a channel of size L x and L y (see 
Fig.[T]). The channel is narrow (L x >> L y ) and at both of 
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its ends there are thermal walls that are kept at tempera- 
tures Ti and T 2 [lj| . After the collision with the wall kept 
at temperature T a disk has its normal component sam- 
pled from the distribution p(v x ) — ^^F^exp(— ^) 
with the sign in the argument of the Heaviside func- 
tion depending on the location of the wall. Its paral- 
lel component is sampled from the Gaussian distribution 

v 2 

p( v v) — "72^r ex P( — 2T )■ ^ n ^ ne vertical y-direction pe- 
riodic boundary conditions are used. To calculate the 
heat conductivity k we use Eq.([T]) with heat flux in the 
x-direction defined as a time average of j = J2i v xi/^ 
where v X i is the x-component of the velocity of the i-th 
particle. 

One can note that without heat reservoirs models 
of this kind are chaotic [HI, Such a feature 

makes this model more realistic than, for example, a 
frequently examined one-dimensional alternating-mass 
hard-disk model [U, [l8| . In the simplest setup of our 
model, known as a 'pen case' model [H, E3 (Fig[lh) 
disks are so large that they cannot exchange their posi- 
tion (Ly/2 > r > Ly/4). Heat conduction was already 
studied in such a case by Deutsch and Narayan [n| . Al- 
though their calculations indicate that a is close to 1/3, 
reported strong finite size effects and relatively small size 
of examined systems (N < 2048) suggest that one has to 
be cautious in the interpretation of these results. When 
r and L y become small and the surface of particles is 
considered as very rough the 'pen case' model becomes 
the random-collision model. Calculations in such a case 
also suggest that a is close to 1/3 fl9l |. 

An efficient way to simulate hard-disk systems is to 
use event-driven molecular dynamics [20l . l2ll . [22j. Per- 
formance of the algorithm considerably increases upon 
implementing heap searching and sectorization and such 
methods have already been applied to a number of prob- 
lems [lH HII • In our model it is sufficient to use sector- 
ization only in x direction. With such a technique we 
examined systems of up to N = 3 • 10 4 hard disks. In the 
'pen-case' version of the model there is no need to intro- 
duce sectorization and simulations are only a little bit 
more demanding than in a one-dimensional alternating- 
mass model 0. 

To allow an exchange of particles we simulated also 
systems with disks of a smaller radius (FigQ})) and in 
such a case the sectorization considerably speeds up sim- 
ulations. To examine the role of momentum conservation 
a fraction c of disks are made immobile. These disks are 
of the same radius and collisions with them conserve en- 
ergy but not momentum. They are placed along the line 
y = Ly/2 (FigQJ) but similar (not presented) results are 
obtained for the random distribution of immobile disks. 
In our simulations the parameters were chosen as follows: 
L x = N, L y = 1.0 - 1.5, r = 0.01 - 0.3, T x = 1, T 2 = 2, 
and c = — 0.1. Initially centers of disks are usually 
uniformly distributed along the line y — L y /2 (for c > 
they are between immobile disks). Their velocities are 
sampled from the Boltzmann distribution at tempera- 
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FIG. 1: In our model disks move between thermal walls kept 
at different temperatures Ti and T2. In vertical direction 
periodic boundary conditions are used. In (a) particles cannot 
pass each other and in (b) such a movement is possible. In 
(c) a fraction of disks is immobile and placed in the middle of 
the system (filled circles). 
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(a) r=0.3,L : 

(b) r=0.1,L;: 
(b) r=0.3, l4=1.5 
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FIG. 2: The size dependence of the heat conductivity re for 
the momentum conserving cases (a) and (b). The solid line 
has a slope corresponding to a — 1/3. For r — 0.3 our results 
for L — 1.0 and 1.5 are nearly identical. 



ture in terp olating linearly between T% at x — and T 2 at 
x = L 25]. Such a system evolves until a stationary state 
is reached and then computations of some time averages 
are made. 

Numerical simulations in the momentum conserving 
case show that k diverges with the number of particles 
N (Fig[2]) and the exponent a is close to the expected 
for one dimensional systems value 1/3. In the limit of 
ideal gas (r — > 0) heat flux is independent on N (note 
that L x = N and increasing the number of particles N 
we also increase the distance they had to travel) This 
explains the slower convergence seen for r = 0.1. 
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(a) r=0.3, L =1.0 

(b) 1=0.3, lL=1.5 

(c) i=0.3, Ly=1.5, c=0.1 
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FIG. 3: The time dependence of the m.s.q. displacement 
{x 2 (t)} for N = 1000 and T 1 = T 2 = 1. 



FIG. 4: The time dependence of the m.s.q. displacement 
{x 2 (t)} for momentum conserving case (a) with different val- 
ues of N and Ti = T 2 = 1.0. 



The divergence of k is an expected feature of momen- 
tum conserving systems @, HJ , and some arguments sug- 
gest that in systems where momentum is not conserved k 
should be finite A simple way to introduce momen- 
tum nonconservation into our system is to place some 
immobile disks. Although we do not present numeri- 
cal data, our simulations confirm that in such a case 
k remains finite in the thermodynamic limit N — > oo. 
We noticed that even a small fraction of immobile disks 
(c = 0.01) is sufficient to remove the divergence of k. It 
means that the behaviour of k in a very sensitive way 
depends on the conservation of momentum. 

Simulations show that when there are no immobile 
disks, k diverges with the same exponent a(~ 1/3) both 
when particles cannot exchange their positions (case (a) 
in Fig fl]) and when they can (case b). If Eq.Q or Eq.Q 
holds, we should observe in both cases the superdiffusive 
behaviour (albeit with different exponents 0). To ex- 
amine diffusive properties we measured the mean square 
displacement (x 2 (t)) over disks that at a certain time 
enter the central part of the system and did not hit a 
wall before the time t (for the examined time scale such 
processes were extremely rare) (26| . To measure (x 2 (t)} 
the system is not subjected to the temperature difference 
(T 1= T 2 =T). 

Numerical results show (FigJ3]) that in the presence of 
immobile disks (x 2 (t)} increases linearly in time ((3 = 1) 
and since k remains finite in this case (a = 0), both © 
and (j4|) are satisfied. When momentum is conserved, two 
different behaviours are observed. In the case (a) disks 
cannot exchange their positions and displacement of par- 
ticles is severely restricted. As a result (x 2 (t)} saturates 
and that is an indication of strong subdiffusion. We will 
see, however, that this is only a finite size effect and in 
the thermodynamic limit a different behaviour emerges 
in this case. In the case (b) disks can exchange their 
positions and asymptotically (x 2 (t)) increases linearly in 
time, as in the momentum nonconserving case. 



Now, let us examine an interesting similarity in the 
behaviour of (x 2 (t)) in the case (a) and (b). Namely, 
initially (x 2 (t)) has some oscillatory behaviour in these 
cases and there is no indication of such a behaviour in 
the momentum nonconserving case. To examine the ori- 
gin of these oscillations we simulated systems of different 
number of disks and at different temperatures. Simula- 
tions show that the time of the first maximum of (x 2 (t)} 
is approximately proportional to N (and thus to L x ) and 
inversely proportional to s/T (note that y/T is propor- 
tional to the typical velocity of disks). Such a behaviour 
indicates that quasioscillations of (x 2 (t)) are related with 
sound-wave effects. Size dependence of these quasioscil- 
lations for the case (a) is shown in Fig[4] and a similar 
behaviour was found for the case (b). Let us also no- 
tice that in case (b) the short-time growth of (x 2 (t)) is 
faster than the long-time growth, although in both cases 
th e growt h is linear. Moreover, the saturation value of 
\J (x 2 (t)} is much smaller than the system length L x . 

Although at the large time scale the behaviour of 
(x 2 (t)} in cases (a) and (b) is much different (FigJ3]), at 
the shorter time scale it shows some similarity. Heat con- 
ductivity in cases (a) and (b) also behaves similarly. It 
is thus tempting to suggest that sound-wave effects, that 
provide a relatively fast transfer of energy but only on 
a short time scale, are related both with divergence of 
k and with quasioscillations of (x 2 (t)). In presence of 
immobile disks such effects disappear, apparently due to 
the dissipation of momentum during collisions with im- 
mobile disks. As a result k remains finite and (x 2 (t)} 
increases monotonously in time. However, more detailed 
studies would be needed to substantiate such a claim. 

Let us also notice that the time scale set by sound-wave 
effects diverges in the limit N — > oo. In that case such 
effects will dominate diffusive behaviour for arbitrarily 
long time. As seen in FigJH in such a limit (x 2 (t)) seems 
to develop longer and longer linear increase. Thus we 
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expect that in the limit N — > oo, the diffusion in both 
cases (a) and (b) (data not shown) is normal, contrary 
to the predictions of (J3]) or J4|. 

That in the case (a) the mean square displacement 
(x 2 (t)} increases linearly in time is perhaps interesting on 
its own. In this case particles cannot exchange their po- 
sitions and that resembles the molecular diffusion e.g., in 
some zeolites [27|. For such, so-called single- file systems, 
the mean square displacement is known to increase as \fi 
and such a slow increase was derived for some stochastic 
lattice gas models [28]. As shown in FigJH continuous 
dynamics and/or momentum conservation considerably 
modify such a behaviour. However, since the thermal 
motion is usually rather fast, the sound-wave time scale 
is quite short and it might be difficult to examine such 
effects experimentally. 



In conclusion, our work shows that in momentum- 
conserving hard disk systems in narrow channels heat 
conductivity k diverges with the exponent a 1/3 and 
sound-wave effects enhance diffusion. As a result, in the 
thermodynamic limit, normal diffusion appears even in 
the 'pen-case' version of our model. When momentum 
is not conserved, k remains finite and no enhancement 
of diffusion was observed. Heat conduction, diffusion, 
and momentum conservation are terms of fundamental 
importance in physics. Their intricate relations even in 
such simple systems like the ones examined in the present 
paper should warrant further study of these problems. 
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